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Periodic  Orbits  near  the  Critical  Inclination  Angle 

by 

Richard  B.  Barrar 

ABSTRACT 

The  existence  of  periodic  orbits  near  the  critical  inclination  angle  for 
the  equations  of  motion  of  a  satellite  of  an  oblate  planet  are  herein 
Investigated.  It  is  shown  that  in  the  immediate  neighborhood  of  the  critical 
angle  these  periodic  orbits  are  librations.  The  period  of  these  orbits  is 
shown  to  be  extremely  long,  i.  e. ,  of  order 

Introduction 

Many  authors  have  treated  the  general  problem  of  the  motion  of  an 
artificial  satellite  in  the  neighborhood  of  the  critical  inclination  angle. 

See  Garfinkel  (1960),  Hagihara  (1961),  Hori  (1960),  Izsak  (1962),  Kozai  (1961), 
for  a  background  in  this  interesting  problem. 

In  the  present  paper  we  restrict  ourselves  to  the  study  of  those  special 
orbits  which  are  periodic. 

The  present  paper  is  an  extension  of  the  author's  previous  paper  on  the 
existence  of  periodic  orbits  away  from  the  critical  inclination  angle  (see 
Barrar  (1963)).  Bie  notation  and  some  of  the  results  of  that  paper  will  be 
used  in  the  present  paper.  For  example,  our  meaning  of  a  periodic  orbit  will 
be  that  of  Barrar  (1963).  (In  Barrar  (1963),  some  results  were  obtained  for 
motion  exactly  at  the  critical  inclination  angle;  the  present  paper  extends 
these  results  to  a  neighborhood  of  the  critical  inclination  angle.) 
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Statement  of  Problem 


We  assume  that  after  a  change  of  variables  (removing  the  short  periodic 

terms)  our  Hamiltonian  can  be  written  in  the  form  (see  Barrar  (1963))^ 
equation  (36): 

F  =  ^^(3^.+  e  cos 


(1) 


where 


(2) 


Ci(ai,a2,Wi,W2)  =  2  A"  cos  (m^^Wj^  +  +  h(aj^,a2)) 


Further  we  assume  (corresponding  to  the  critical  angle)  that  there  exists  a 
point  such  that: 

(3)  nj^  =  dCj^(apa°)/da^  ^0  n^  =  ^j^(a°,a2)/5a2  =  0 

but 

A  =  b^C^(,a°^,a°^ /'da^ba^  ^0  D  =  S^Cj^(a°,a°)  ^aj^oa^^  ^  0 

(4) 

B  =  a^Cj(a°,ap/aa^da2  =  0 

Since  n^  vanishes,  the  method  developed  by  Poincare  (1893  §134)  (i. e. , 


the  von-Zeipel  method)  for  non-critical  angles  no  longer  applies.  However, 
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Poincare  (1893  §206)  develops  another  method  that  is  applicable  to  the  Hamiltonian  (1) 
together  with  (3)  and  (4).  This  method  converts  the  Hamiltonian  (1)  Into  another 
Hamiltonian  whose  form  Is  the  same  as  that  of  the  original  Hamiltonian  but 
with  the  corresponding  n|^  and  n^  not  vanishing.  Poincare  calls  this  procedure 
the  method  of  Bohlln.  We  give  a  presentation  of  this  method  below^  applicable 
to  the  Hamiltonian  (1)^  (2)^  (3)  and  (4). 

The  Polncare-Bohlln  Transformation 

As  motivation  for  the  Polncare-Bohlln  transformation  let  us  first  consider 
the  Hamiltonian: 

^*^®l^®2^"l'V  “  ‘^l^^l^V  S2^®1'*2^  ^''2 

(5) 

+  C23(a^,a2,WpW2)] 

together  with  the  conditions  (3)  and  (4).  We  Introduce  S  =  S  +  S,  + 

o  1 

€82  with  *  a°  Wj^  +  a°  and  wish  to  solve  the  equation: 

(6)  F*(dS/^Wj^,dS/dw2,Wj^,W2)  -  +  e  D2  +  0(e^/^) 

1/2 

If  like  powers  of  €  '  are  equated  in  this  equation,  one  obtains  the  series  of 
equations: 

O)  c^(a°,a°)  -  D^(a;,a°) 

nj  =  Dj^ 


(8) 
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(9) 


+  °2  ■  • 


It  follows  from  this  last  equation  that  in  order  for  $2  not  to  have  any  secular 
terms  one  must  have: 


(10)  +  D(aSj^/Swp^  “  °2  ■  Sl^®1^^2^  ■  ^22^^V^2^  ^^^2 

*  1/2 

Now  if  we  solve  only  (7),  (8)  and  (10)  and  set  S  =  +  e  '  we  obtain 


(11) 


F*(as*/Sw^^  ^S*/^2^  ''l^''2^  “  °o^®P®2^  ■*■  °1 

+  6(D2-C23(a°,a°,w^,W2))+0(e3/2) 


o  o 


In  the  above  constants.  If  in  terms  of  new  variables  we  can 

make  Dj^  =  Dj^(xp,  D2  =  D2(xpxp  with  dD2/dx^  ^  0,  then  in  terms  of  these  new 
variables  we  are  no  longer  faced  with  the  critical  inclination  angle  difficulties. 
We  do  this  more  precisely  below.  (The  above  is  a  slight  modification  of  the 
method  used  in  Poincare  (1893)  §206.  Since  Poincare  treats  a  Hamiltonian  of  the 
form  (5),  but  we  must  treat  one  of  the  form  (1),  this  modification  is  essential.) 

We  now  transfer  from  the  variables  aj^,a2, Wj^,W2  to  new  variables  xpx^y^y^ 
by  means  of  a  generating  function 

(12)  T(a°,a2^xpxpWj^,W2)  =  +  a2'^2  "*■ 


with  satisfying: 
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(13)  »  x| 

(14a)  A(dTj^/dw2)^  +  D(xp^  “  '  ^21^®1'®2^  "  ^22^^1^®2^  ^'^2 

Further  set: 

(14b)  Sl^®l^®2^  ■*■  ^22^®1^^2^  2w2  =  A  t 

(14c)  =  (D^  -  D(xp^)/A 

Then, 

(15)  Tj^  *  xj^w^  +  J  2  (5^2 

[Thus,  in  (11)  *  n^xp  »  A(x^  +  D(xp^).  ] 

The  canonical  variables  sre  determined  from  the  equations: 


(16) 

Sj,  =  St/ow^  =  Sr/dxj 

or 

(17a) 

Sj^  =  a°  +  x^ 

(17b) 

y{  * 

(17c) 

a2  -  a°  +  e^/^(x^-t)^/^ 

(17d) 

y^  =(€^/V2)r(x^-,^)'^/2 
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[For  the  physical  meaning  of  this  transformation  see  Poincare  (1893)  §199^  §200,  and 
Hagihara  (1961),  discussion  of  Figure  1.  ] 

Because  of  the  form  of  t  (see  (14b)),  integral  (17d)  is  an  elliptic  integral, 
and  sin  w^,  cos  are  doubly  periodic  in  y^.  If  is  greater  than  the  maximum 
of  \lr,  the  real  period  is 

2n 

(18a)  €^/^P(xp  =  (xpt) 

However,  if  xpt  vanishes  for  =  a,  ^  and  remains  positive  for 

cc  <  <  P,  then  the  real  period  is 

P 

(18b)  €’'^^P(xi)  =  (e^^Vz)  f 

^  '^OL  ^ 

In  the  case  (18a)  increases  by  211  when  y^  increases  by  a  period.  In  the  case 
(18b),  that  is  in  the  case  of  libration,  returns  to  its  primitive  value  when 
y^  increases  by  a  period.  Hence,  the  variables  x^x^y^y^  are  uniformizing 
variables  for  the  original  variables  a^^, a2,Wj^,W2. 

Several  interesting  examples  arise  in  the  motion  of  an  artificial  satellite 
near  the  critical  inclination  angle  as  noted  by  Hori  (1960).  First,  when  only 
the  second  harmonic  is  considered  then  only  (18a)  applies.  Secondly,  when  the 
Vinti  (1959)  potential  is  considered  the  terms  involving  cosines  and  sines  of 
in  ^  vanish  and  hence  (17d)  is  essentially  a  linear  transformation,  with 
(18a)  applying.  Finally  when  both  the  second  and  fourth  harmonics  of  the  earth's 
potential  are  considered,  then  for  sufficiently  small  values  of  xp  the  case 
of  libration  (18b)  applies. 
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If  one  makes  the  change  of  variable  (17)  in  the  Hamiltonian  (1),  and  takes 
(11)  into  account,  he  finds: 


*^22  (aj^,a2)  cos  2w^ 


(19) 


C  (ap  +  e’-Z^F* 
o'  1^ 


with 


F*  =  jy'ji  ^  g  +  e^Z^A(x^  +  D(xp^) 

I 

+  €^/^C23(a°,a°,w^,W2)  +  O(e^) 


where 


C^^^ap  =  d^C^(a°)/(dapj 


Further  if  one  introduces  the  change  of  variable 


(21) 


y'  =  z^e 


1/2 


and  since  C^(aj^)  In  (19)  is  a  constant  the  equations  of  motion  become: 


^^^The  present  method  of  dealing  first  with  (5)  and  then  with  (1),  in  making 
change  of  variable,  was  used  to  avoid  the  difficulty  mentioned  in  the  introduction 
of  Izsak  (1962)  with  terms  of  order 
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(22)  dx^/dt  =  dz^/dt  = 

Although  F  is  periodic  with  respect  to  the  period  with  respect  to 

z^  is  2IIj  but  with  respect  to  z^  it  is  P(xp  as  defined  in  (18).  To  make  the 
period  211  with  respect  to  both  angle  variables  we  introduce  the  last  canonical 
variable  change; 

(23)  u^  =J  P(xpdx^/2n  =  2nz2/P(xp 

The  equations  remain  canonical; 

dx|/dt  =  hF*/bz^  du^/dt  =  b?*/bv^ 

(24) 

dZj^/dt  =  -  dF*/dx|  dv^/dt  =  -  dF*/du2  . 

F  is  periodic  with  period  211  with  respect  to  original 

variables  aj^,a2,w^,W2  can  be  expressed  in  terms  of  x^,U2^ Zj^, and  are  periodic 
with  period  2II  with  respect  to  Zj^^v^. 

Periodic  Solutions 

By  means  of  the  Poincare-Bohlin  transformation  we  have  changed  the  Hamiltonian 


(1)  to  one  of  the  form; 


28  January  1963 


9 


SP- 106 7/001/00 


where 

£^2  =  2  A'(x|,U2)  cos  +  hCx^^u^)) 

ra^/0 

Ej  =  E  A"(x|,U2)  cos 

with 

d  Eo(xp/dxj[  ^  0  dE2^(x^,U2)/du2  /  0  . 

This  Hamiltonian  (25)  is  slightly  different  than  the  Hamiltonian  treated  by 
means  of  the  Poincare — von-Zeipel  transformation,  in  Barrar  (1963),  in  that  the 
variable  U2  does  not  occur  until  the  fourth  term.  Nevertheless  the  Poincare-- 
von-Zeipel  transformation  is  applicable  to  the  Hamiltonian  (25)  with  some 

1/2 

obvious  modifications  such  as  expanding  in  powers  of  €  '  not  e.  Thus  from 
the  results  of  Barrar  (1963),  we  conclude  that  (25)  has  periodic  solutions  for 
arbitrary  initial  values  of  and  for  values  of  satisfying  the  Poincare 

criteria  (19b)  of  Barrar  (1963).  Since  the  original  variables 
periodic  in  x|,U2,z^,V2  this  also  insures  periodic  solutions  in  terms  of  the 
original  variables. 

Note  that  the  period  of  the  resulting  motion  will  be  very  long  since  the 

3/2 

frequency  with  respect  to  V2  is  0(e  '  )  of  the  frequency  with  respect  to  Zy 
Conclusion 

We  have  studied  periodic  orbits  in  the  immediate  neighborhood  of  the  critical 
inclination  angle  for  the  motion  of  a  satellite  of  an  oblate  planet.  By  means  of 
the  Poincare -Boh 1 in  transformation  we  have  been  able  to  introduce  new  uniformizing 
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variables^  such  that  the  equations  of  motion  In  terms  of  the  new  variables  have 
essentially  the  same  form  as  the  equations  of  motion  In  terms  of  the  old  variables 
but  away  from  the  critical  angle.  From  the  results  of  Barrar  (1963)  we  are  then 
able  to  conclude  the  existence  of  periodic  orbits  near  the  critical  Inclination 
angle.  Very  near  the  critical  angle  these  periodic  orbits  are  llbratlons. 
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